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Abstract. - We study the fluctuations of the two-time dependent global roughness of finite size 
elastic lines in a quenched random environment. We propose a scaling form for the roughness 
distribution function that accounts for the two-time, temperature, and size dependence. At 
high temperature and in the final stationary regime before saturation the fluctuations are 
as the ones of the Edwards- Wilkinson interface evolving from typical initial conditions. We 
analyze the variation of the scaling function within the aging regime and with the distance 
from saturation. We speculate on the relevance of our results to describe the fluctuations of 
other non-equilibrium systems such as models at criticality. 



The study of dynamic fluctuations may serve to grasp the nature of non-equihbrium phe- 
nomena. Building upon previous analysis of critical static and dynamic fluctuations [1,2], 
Racz [3] proposed to use the scaling functions characterizing the fluctuations of macroscopic 
- global - observables in non- equilibrium steady states as a classification tool. By introducing 
'universality classes' in this way, one could then use them to uncover symmetries and dynamic 
mechanisms in experimental systems. 

Another class of non-equilibrium phenomena, including important problems such as glassi- 
ness and coarsening, is characterized by a slow relaxation with loss of stationarity [4,5]. It is 
by now becoming clear that the study of fluctuations [6] is necessary to understand the mech- 
anism for the slowing down and the aging non- equilibrium relaxation in these cases too [7,8]. 

An elastic line under the effect of quenched disorder is a relatively simple system with many 
aspects of glassiness due to the competition between elasticity and disorder. Its physical real- 
izations are manifold, including interfaces in random ferromagnets [9], crack propagation [10], 
and vortex lines in high- Tc dirty superconductors [11]. The global noise and disorder averaged 
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displacement and linear response age with a multiplicative temporal scaling [12, 13] similar to 
what is found at criticality [14] and in Sinai diffusion [15]. In this Letter we study the line's 
fluctuating relaxation by analyzing the probability distribution functions (pdfs) of the two- 
time dependent global roughness and we discuss it in comparison to what has been found in 
other aging systems [6-8] and the Edwards-Wilkinson (EW) interface [2]. 

We study a lattice string of length L directed along the y direction of a rectangular square 
lattice of transverse size M ^ L"^/^ [M — 10^, L = 500) ensuring the existence of many nearly 
equivalent ground states. The line segments, Xy {y = 1, - ■ ■ ,L), obey the restricted solid-on- 
solid (SOS) rule \xy—Xy-\\ = 0, 1. A quenched random potential V taking independent values 
on each lattice site is drawn from a uniform distribution in [—1, 1]; we use 10^ — lO'^ realization 
depending on L. The initial configuration {t = 0) is drawn from the equilibrium distribution 
at the high T = 5. At each microscopic time step we attempt a move of a randomly chosen 
segment to one of its neighbours restricted by the SOS condition and we accept it with the 
heat-bath rule. One Monte Carlo (MS) step is defined as L update attempts. In the following 
we use adimcnsional time, space and energy scales [12]. 

The glassy phenomenon appears in this model as a dynamic crossover [12,13]. For all 
observation times, toi,g. that are longer than a size and temperature dependent equilibration 
time, teq, the dynamics is stationary. Instead, for tobs < teq the relaxation occurs out of 
equilibrium as demonstrated by two-time correlations that age. The latter are measured as 
follows. After equilibration at high T the system is quenched to low T and time is set to 
zero. The line relaxes until a waiting-time, t^,, when the quantities of interest are recorded 
and later compared to their values at a subsequent time t. 

In this Letter we focus on the two-time global roughness 

W\t,t^) = L-^T.y^A53:y{t) - 5Xy{t^)\^ , (1) 

with 5xy{t) = Xy{t) — x{t) the distance from the center of mass, x{t) = L^^ S^=i ^yi^)- Before 
any averaging, vcP' fluctuates when changing the disorder and thermal noise realizations. 

At high T, neglecting a short transient, the disorder and thermal averaged roughness, {w^), 
is stationary and crosses over from growth to saturation at tx ~ [16]: 

{w'^{t,t^)) r^L^ f{At/tx) , At = t-t^, (2) 

where the scaling function obeys f{x) ~ x^ for x <^ 1 and f{x) ~ 1 for a; ^ 1, with (3 
and z = the roughness, growth and dynamic exponents, respectively. For the EW line 
in 1 + 1 dimensions C = 1, /? = 1/2 and z = 2. In the presence of disorder, C is expected to 
take a 'thermal' value, Ct/n for L < L^iT), and a larger 'disorder' dominated value, Qis, for 
L > Lc{T), both being T-independent [16]. 

For not too short L and sufficiently low T, teq goes beyond the numerically accessible 
time- window and (w^) depends on tw explicitly. Before saturation one expects 

{w^t, t^)) ~ 6{t^) {w^[i{t)ll{t^)]) = f{t^) Tm/i{t^)] (3) 

with l{t) a growing length (dimensions are restored by prefactors that we omit). This form 
approaches a stationary regime when -C At <C tx if J^{x) x"^ iov x ^ 1, and complete 
saturation at when i{t) —* L. We found data collapse using i{t) ~ t"/^ with a small 
exponent [17] and r[l{t) / t{t^)\ = g{At/ty,) with 



g{x) ~ x'^A lO^s(^) ^ x'^A 10 



(4) 
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Fig. 1 - (Colour on-line.) Panel (a): thermal and disorder averaged square roughness as a function 
of time-difference for a line with L — 500. The top curve is for T = 5; the dotted lines indicate the 
sub-diffusive regime, saturation and the characteristic time tx- The high-T dynamics is stationary and 
conforms with the scaling in The bottom curves are for T = 0.2, and = 10^ 10^ fO", 10^ lO'' 
MCs from top to bottom. At low T, stationarity is lost after an initial At regime. Saturation is out 
of the At window for all low T curves. Panel (b): sketch of the time-difference dependence of the 
averaged roughness for several waiting-times. The solid lines have been drawn using the function Q in 
eq. (2J and a crossover to saturation at (it)^)oo. The dashed lines represent the asymptotic power laws 
ci,2(r)At°. The horizontal dotted line indicates a path of constant (w^) / (uP')ac,- The intersection of 
the curves with the inclined dotted line correspond to constant u\, see eq. O. Panel (c): scaling of 
the data in panel (a) using eqs. Q-Q; the solid line is 5 in @ with a = 0.145, A = 0.28, B = 0.35, 
C = 0.56 and D = 5.01. Panel (d): L-dependence of the saturation value, (w^)oo, at high and low T. 

in the limits a; 1 and x ^1, respectively, and g{x) = tanh[Clogjo(^/-'^)] which does not 
have a special significance but serves to select the working times At and i^, below. Note that 
on the two asymptotes stationarity is recovered: 

(w^) ^ ci,2(T) At" , with a < Pew = 0.5 , (5) 

and different proportionality constants ci^2 = A lO^'^. a is then a generalization of the growth 
exponent, f3 and a < Pew reflects that disorder slows down the dynamics. In Fig. ^ we show 
{up') at high and low T: panels (a) and (c) display numerical data and panel (b) presents a 
sketch of the size and time-dependence of {vP) at low T. Details are given in the caption. 

Even though (w^) weakly depends on the polymer size before saturation, the saturation 
value strongly depends on L and it turns out to be important to describe the fluctuations. 
Transfer matrix calculations show that the static roughness Wy = {[xy — xqY) with xq fixed 
crosses over from Wy ~ y to Wy ~ y'^^^ at a length y* '--^ 5 at T = 0.2 [12]. One might expect 
a similar crossover for the asymptotic global roughness ^ though the crossover length y* 
should not necessarily be the same. With MC simulations we can only determine (u'^)oo = 
limAt>t;j (w^) for rather short polymer lengths {L < 32) both at high and low temperature. 
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Fig. 2 - (Colour on-line.) Bare square roughness fluctuations at T = 5 [(a) and (c)] and T = 0.2 
[(b) and (d)] for a line with L — 500. The time-differences are given in the key. The waiting-times 
are tm = 10^, 10^, 10^, 10^ MCs [from right to left in panels (b) and (d)]. In panel (a) one curve 
(dashed green) for L = 100 and the same parameters is included to highlight the L-dependence of 
the fluctuations. 

Figure ^d) shows its i-dependence at T = 5 and T = 0.2, two values that we shall use in 
the rest of the Letter. The dependence is linear with no sign of the crossover to the disorder- 
dominated regime up to the maximum length, L = 50. 

We now turn to the study of fluctuations. In Fig. [3 we show the pdf of square roughness 
fluctuations, p{w'^), for a line with L = 500 at T = 5 [left, panels (a) and (c)] and T = 0.2 
[right, panels (b) and (d)] for several values of and At given in the keys. At T = 5 there is 
no tw dependence. The pdfs evolve with At and they basically take the EW form [2]. At low 
T a clear dependence appears at sufficiently (but not too) long At's. 

In order to understand p{w'^) wc propose a scaling form that we later put to the test 
numerically. In the low T regime, for times such that equilibration has not been reached 
{tuj <§; t^g) p{w^) depends on the (already adimensional) times t and t^, the system size L 
and temperature T. With a simple change of variables we measure w'^ in units of its average 
(ui^): p{w^;t,t^, L,T) — {w^)^^p{w'^ / {w'^);t,t^, L,T). Using eq. ©, and the monotonicity 
of £ and Q, the dependence on t and can be replaced by a dependence on (w^) and (w^). 
L can be traded for the T-dependent saturation value (i«^)oo- Now, our scaling hypothesis is 
that the dependence on (w^) occurs only in comparison with the saturation value: 

{w^)piw') = $ (^^'), {w^)/{w^)oo,T) . (6) 

This form holds exactly at high temperatures and for the EW interface. In these cases dis- 
order is irrelevant or even absent and the f^-dependence, and hence the (■£'^)-dependence, 
disappears. In addition, all T dependence enters through (w^)oo and one has {'uP')p{nP') = 
^Ew{wy{w^),{w^)/{w^)^) [2]. 
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and 1 when {uP') moves from the first to the second asymptotic power laws before saturation, 
see the dashed lines in Fig.^b) and eq. IS}- We call V2 = (w^)/ (u'^)oo the second parameter 
in ©. Then, by choosing groups of (L, Af, ty,) we vary vi and 1^2- In the following we test 
the following restatement of our scaling hypothesis 



In Fig. I^Ja) we plot data for i/i = 0.389 and different i'2- We use (w^)oo ^ L and we plot 
six sets of data corresponding to two values of 1^2 = {w'^)/L (the T-dependent prefactor is 
irrelevant since we are working at constant T) . The two sets of data with the same V2 collapse 
rather well. The dependence on V2 is similar to the one found for the EW line: for small V2 
the scaling function has an approximately log-normal form, it is positively skewed and it gets 
wider with increasing V2. In panel (b) we use instead (w^)oo L^^^ and the scaling variable 
V2 — {'uP') / L'^/^ . The collapse is now lost in accordance with our hypothesis and the fact that 
{w'^)oo ~ L for the lengths used. (It would of course be desirable to have an independent 
determination of (w^}oo to confirm our conclusions.) 

In Fig. |3{c) we analyze the evolution of $ with vi. $ resembles the high T result with 
the same distance from saturation (black solid line) when vi ^ \ and {np) ~ C2(r)AP. It 
progressively deviates from the high T form for decreasing vi . Decreasing ui while keeping V2 
fixed has a similar effect as increasing V2 while keeping ui fixed: $ gets closer to the asymptotic 
disorder-free form, limAt>t^ '^ew ■ Moreover, at fixed and small V2 the j/i-dependent pdfs can 
also be quite well described with a log-normal function with two i^i-dependent parameters. 

The vi dependence can be rationalised from the point of view of the dynamics in a free- 
energy landscape. For each there is a sufficiently long At such that {vo^) ^ C2(T)At" and 
the waiting-time dependence disappears (see Fig. QJ. After these very long At's one allows 
the system to get out of the traps occupied at t^,, sub-diffuse with a < (3ew — \ but fluctuate 
similarly to the high T and Gaussian cases once the good scaling variables are chosen, as shown 
by the fact that p{'uP') gets very close to the continuous curve in Fig.lS^c). In the opposite limit 
in which {uP) has not deviated much from the first asymptote the system is quite trapped 
and has not had enough time-difference. At, to explore a relevant part of the free-energy 
landscape. The fluctuations are then very different from the high-T and disorder-free ones. 

In the inset to Fig. Efa) we show data from the Langevin dynamics of a continuous model 
of an elastic line in a Zd disordered environment [13, 20] that confirms the scaling form and 
the generic trend. We shall present more details on the functional form of $ in the lattice and 
continuous models as well as on the T-dependence in a longer publication [23]. 

Summarizing, wc studied the effect of disorder on the low-T averaged and fluctuating two- 
time roughness of elastic lines. On the one hand disorder introduces aging at sufficiently low 
temperatures. Thus, a time-delay that increases with the waiting-time is necessary to enter 
a stationary regime next to saturation. On the other hand, it slows down the evolution and 
the averaged roughness in this last dynamic regime undergoes sub-diffusion with a smaller 
exponent than in the EW case. Well before saturation the two-time roughness fluctuations 
can be represented by a scaling function, that depends on where in the aging regime the 
dynamics take place (parametrized by {uP) or vi) and the distance from saturation (given 
by V2 = {up) / {up) 00) ■ At high T, and at low T but close to the sub-diffusive regime before 
saturation, disorder is quite irrelevant and <f> approaches the EW form ^ew depending only 
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Fig. 3 - (Colour on-line.) Panel (a): check of the scaling Q using (w^)oo oc L (C = !)• The gray 
filled squares (L = 100, — 100), pink open squares {L = 200, fu, = 1.19 10*) and cyan open circles 
{L = 300, t^o — 1.95 10^) data are at the same distance from saturation {v2 = {'uP)/L — 8.8 10~^), 
and have the same 'aging' roughness {u\_ — 0.389, At/tw = 10). The blue -k (L — 200, = lO'^), 
green x (L = 400, t» = 1.19 10^) and red + {L = 600, tn, = 1.95 10^) data are for U2 = 4.4 10"^ and 
the same as the other group. Inset: the pdf for the continuous model in [13]. Panel (b): test of the 
scaling @ using (w^)oo oc L*^^ {( = 4/3) and U2 = (w^)/L'^^^ = 1.3 10"^. Panel (c): Dependence on 
vi and comparison with the high-T behaviour (black solid line). One gets deeper into the aging regime 
in the order blue ★ (At = 2.25 10^^» = 10^ i^i = 0.09), green x {At = 1.50 10^^l„ = 10^ i/i = 0.44), 
red + {At= 1.50 10^^„ = W°,ui = 0.72). U2 = = 2 10"^ 



on 1/2- At intermediate time-scales, where the waiting-time dependence is explicit, $ deviates 
from the free result and becomes broader and less symmetric. The dependence on the distance 
to saturation at fixed 'aging' parameter z/i is similar to the one in the EW Une. 

We conjecture that this kind of scahng appears also in line problems in different universality 
classes, e.g. KPZ [21], and/or for other quantities such as the maximum height with respect 
to the mean [22]. It would also be interesting to study these features in the out of equilibrium 
dynamics of critical systems [14] and random walkers in random environments [15] with a 
multiplicative scaling of two-time correlations and responses. 

This study extends the proposal in [7] for the fluctuations in aging systems with additive 
scaling - correlations that approach a well-defined and finite plateau to later further decay to 
zero - to problems with diffusive aging and a strong dependence on the system size. In the 
former cases one could explain the scaling and generic form of the pdfs as a consequence of the 
emergence of a symmetry, time reparametrization invariance, in the long time dynamics. The 
presence of the symmetry allowed one to propose an effective action for the most important 
fluctuations, and derive from it the pdfs of local correlations and responses that have been 
partially confronted to experimental measurements in colloidal suspensions [6] and simulations 
of kinetically constrained models and finite dimensional spin-glass models [7] . The study of 
the symmetry properties of the long-times dynamics of elastic lines in random environments 
and further identification of an effective sigma-model-type action for the fluctuations remains 
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to be done. It should be possible to use it to compute 
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